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Summary. 

We present results for the 1 dimensional stochastically forced Burgers 
equation when the spatial range of the forcing varies. As the range of forc- 
ing moves from small scales to large scales, the system goes from a chaotic, 
structureless state to a structured state dominated by shocks. This transi- 
tion takes place through an intermediate region where the system exhibits 
rich multifractal behavior. This is mainly the region of interest to us. We 
only mention in passing the hydrodynamic limit of forcing confined to large 
scales, where much work has taken place since that of Polyakov^]. 

In order to make the general framework clear, we give an introduction to 
aspects of isotropic, homogeneous turbulence, a description of Kolmogorov 
scaling, and, with the help of a simple model, an introduction to the lan- 
guage of multifractality which is used to discuss intermittency corrections 
to scaling. 

We continue with a general discussion of the Burgers equation and forc- 
ing, and some aspects of three dimensional turbulence where - because of 
the mathematical analogy between equations derived from the Navier-Stokes 
and Burgers equations - one can gain insight from the study of the simpler 
stochastic Burgers equation. These aspects concern the connection of dissi- 
pation rate intermittency exponents with those characterizing the structure 
functions of the velocity field, and the dynamical behavior, characterized 
by different time constants, of velocity structure functions. We also show 
how the exponents characterizing the multifractal behavior of velocity struc- 
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ture functions in the above mentioned transition region can effectively be 
calculated in the case of the stochastic Burgers equation. 
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I. Introduction. 

We study some aspects of statistically stationary, homogeneous and 
isotropic fully developed turbulence. This is the typical framework in which 
such studies are done. The quantities of interest are the equal time spatial 
correlations of the velocity field u(x,t), the so-called structure functions. 
The longitudinal structure functions, which are the ones usually discussed, 
are defined by 

S P ( r ) =< [(u(x + f,t) -u(x,t)).n] p > (1) 

where n is the unit vector in the direction r. Some components of the velocity 
field difference can be projected onto the direction transverse to n, and thus 
there are other correlations of p — th order, which involve longitudinal and 
an (even) number of transverse projections. 
The velocity satisfies the incompressible Navier-Stokes equation 

d t u + u.Vu = -Vp + vAu + f (2) 

with 

V.u = (3) 

Here p is the pressure divided by the constant mass density, v the kinematic 
viscosity. We have added / = f(x,t), an external stochastic force which 
acts on large scales, and maintains a turbulent steady state. The average in 
(1) then includes as well an average over time. 

In the usual picture of turbulence (see I.I.), when the distance r = \r\ in (1) 
is small compared to large scales L of the order of the system size, and large 
compared to the scales where dissipation takes place, the structure functions 
are expected to behave as 

S p (r) ~ (r/L)^ (4) 

An important aspect of solving the problem of statistical isotropic, homoge- 
neous turbulence is deriving the values of the exponents C p in (4) from the 
Navier-Stokes equation. This has not been done except for £3, the value of 
which is fixed by the Von Karman-Howarth relation It turns out how- 
ever that the experimentally measured C p 's (up to p = 10 or so) are not 
too different from their scaling values as they arise in the picture of fully 
developed turbulence proposed by Kolmogorov. This is the reason a large 
number of phenomenological models exist, which by breaking scale invari- 
ance slightly, give improved fits to the data. The usual language in which 
to express deviations from scaling is that of multiscaling or multifractality. 
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We will therefore discuss first in this introductory section Kolmogorov 
scaling, then a simple model, which allows one to introduce non-scaling 
elements, and provides a simple introduction to the language of multiscaling 
which we present next. A general reference for these subjects is the book of 



In the second section we discuss the stochastic Burgers equation, its 
shock structure and the associated extreme multifractality, and its behav- 
ior when the spatial range of the random forcing varies from small to large 
scales. In section III we take up the point about statistical aspects of the 
stochastic Burgers equation and their connection with three dimensional, 
forced, isotropic and homogeneous turbulence. First we show how the prob- 
lem of multifractality can be solved for the stochastic Burgers equation. 
Then we discuss the relation between intermittency in the energy dissipa- 
tion to intermittency in the velocity field, and end up by making a number 
of observations concerning the dynamical behavior of structure functions. 
General remarks about intermittency in fully developed turbulence and for 
the stochastic Burgers equation are made in section IV. 

This report is based on a number of results or points made in references [Q, 



1.1. Kolmogorov scaling. 

The picture is that of an energy cascade from the large scale L where 
the energy is put into the system, to the dissipation scale 5 where it is 
dissipated. On intermediate scales 5 <C r <C L, which make up the so- 
called inertial range, the only quantity which matters is e, the mean energy 
dissipation rate per unit mass, considered to be independent of scale, e has 
the dimension of velocity squared divided by time, or velocity cubed divided 
by distance. 

The dissipation scale 5 can only depend on v and e, and thus for dimen- 
sional reasons 5 ~ (^ 3 /e)3 ~ (l/Re)^L, where after replacing e in terms of 
a characteristic velocity U and the large scale L, we are able to introduce 
the Reynolds number Re = UL/v. In the limit of small viscosity or large 
Reynolds number there is thus a definite separation of scales between 8 and 



In the inertial region, dimensions are determined by e alone, and therefore 
one predicts on dimensional grounds, that S p (r) which has the dimension of 
velocity to the p — th power behaves as 



Frisch||]. 
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(5) 
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This is Kolomogorov scaling. The scaling values of the exponents in (4) are 
then 

C P =p/3 (6) 

This gives £2 = 2/3, which by Fourier transform is equivalent to the exper- 
imentally observed —5/3 behavior of the energy spectrum, namely E(k) = 
k 2 < u(k).u(—k) >~ esk~3. One also obtains £3 = 1, which is the value 
fixed by the Von Karman-Howarth relation. The other general result Q is 
that C p is a convex function of p. Measurements of the structure functions 
showQ however that Kolmogorov scaling does not hold: the measured £ p 's 
for p > 3 lie below the scaling values. For instance = 1-80 ± 0.05 rather 
than the scaling value of 2, obtained from (6). This effect is called intermit- 
tency or multifractality, and can be related heuristically to the non-space 
filling property of the eddies which make up the energy cascade, and there- 
fore to their fractal dimension. A simple model will serve to illustrate these 
points. 

1.2. A simple model. 

Among models which describe the energy cascade, the so-called /?-model 
| |10| | is instructive. Imagine, as the energy cascades down to smaller scales 
from the large scale L, that at scales r = a n L in the inertial range, the eddies 
at this scale, which themselves have a typical size of r, occupy only a fraction 
j3 of the available space, such that p r = (3 n , where p r can be interpreted as 
the probability of finding an eddy of size r at scale r. Eliminating the 
"generation" number n between the expressions for r and p r , on finds 

p r = {r/Lf~ D (7) 

where 3 — D = ln/3/ln a. If the eddies are space filling, then (3=1, and 
therefore D = 3. The value of 3 corresponds to the fact that we pretend our 
discussion is for eddies in 3 dimensions. The argument itself is clearly inde- 
pendent of space dimension. One now interprets D as the fractal dimension 
of the space on which the eddies exist, assuming that D is smaller than 3. 
What are the structure functions in this model? 

The typical energy of an eddy of size r is E r ~ 5v r 2 p r , and therefore 
the average energy dissipation rate (per unit mass) at scale r, with a typical 
time scale t r = r/5v r , is 

e r ~ 5 -^{r/Lf- D -' (8) 
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Here 5v r is the velocity variation across the eddy. The value of e r is indepen- 
dent of r if homogeneity holds (existence of an inertial scale) , and therefore 
one has for the velocity 

5v r ~ (eL) 1 / 3 ( r /L)i-( 3 - D )/ 3 (9) 

from which follows for the structure function 

5 p (r) =< 5v r p >= SvP Pr ~ ( e L) p / 3 (r/L) p / 3+ ( 3 ~ D )( 1 - p / 3 ) (10) 

One thus finds for the exponents £ p of the structure functions, a convex 
function of p, namely 

C p =p/3 + (3- J D)(l-p/3) (11) 

which satisfies the condition (Von Karman-Howarth relation) £3 = 1. The 
scaling violating part in £ p is given by (3 — D)(l — p/3). For instance = 
2 — (3 — D), which, by comparison with the experimental result £ 6 = 2 — 0.2, 
leads to a fractal dimension D = 2.8. Note that the velocity variation at r 
(5v r ~ r h ) is itself characterized by an exponent h = 1/3 — (3 — D)/3. For 
D = 3, when the eddies fill all space at any inertial scale, one has the scaling 
(Kolmogorov) result h = 1/3 and ( p = p/3. 

In the simple model we have considered, the structure functions and 
the variations of the velocity field are characterized by a single h and D. 
However here, as opposed to the Kolmogorov scaling behavior, the eddies 
are not space filling, but are characterized by a fractal dimension D. 

Simple fractal models such as the one we have described are not believed 
to give the whole picture required to describe fully developed turbulence. 
Experimental data suggest that £ p depends non linearly on p in contrast to 
equation (11). It is believed(3j that one needs to consider a more general 
picture, with a range of possible /i's and of corresponding fractal dimen- 
sions D{h) (see section IV.). This picture, or the language in which it is 
formulated, is that of multifractality, which we discuss next. 

1.3. The language of multifractality. 

Assume now that h can take values in an interval (h m i n ,h max ), and 
that to each h there corresponds a set in three dimensional space of fractal 
dimension D(h), in such a way that across any distance r ( r belongs to the 
inertial range) in the vicinity of that set, one has 

5v r ~ (r/L) h (12) 



6 



and 



Pr ~ {r/Lf- D ^ 



(13) 



where p r is the probability for being within a distance of the set of fractal 
dimension D(h), and 5v r is the velocity variation. As a consequence one has 
the following expression for S p (r) for a given set with scaling dimension h 

S p (r) ~< 5v r p >~ (r/L) ph+3 ~ D ^ (14) 

All /i can contribute to the right-hand side, but since r/L -C 1, the dominant 
exponent ( p is given by 

C„ = mm(ph + 3 - D(h)) (15) 

h 

This exponent Q p is the dominant one in the expression of the structure 
factors (cf. equation (4)). 
Remarks: 

- the scaling result corresponds to h = 1/3 and .0(1/3) = 3. 

- the argument is the same in 1 or 2 dimensions with the replacement of the 
number 3 in 3 — D(h) by respectively 1 and 2. 

- the quantity 3 — D(h) is positive or zero, since D(h) cannot exceed the 
dimension of the embedding space. It is generally assumed that h m i n > 0. 
In the case of the Burgers equation where exponents can be calculated, we 
find (cf. section III.l.) that the h's corresponding to higher order structure 
functions reach the value when the stochastic forcing has moved to suf- 
ficiently large scales, and stay at the value when the scale of the forcing 
increases further. 



II. The stochastic Burgers equation. 



This is a 1 dimensional version of the Navier-Stokes equation, a version 
without incompressibility and pressure, which describes the evolution of the 
compressible field u(x,t), by 

d t u + u— = u—r + / (16) 
ox ox z 

where / = f(x,t) is a stochastic forcing. 

We will discuss later the forcing and its influence on the dynamics of 
the field. For the moment, we will ignore it, and summarize some results 
concerning the plain Burgers equation[ll|. 
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II. 1. Shock structure and extreme multifractality. 

If one starts from an initial sinusoidal velocity profile of large wavelength, 
then under the influence of the nonlinear term in the equation, the sinusoid 
will for sufficiently small viscosity, steepen into a series of shocks. After 
some time the shocks will fade away, their energy being dissipated by the 
viscous term. This viscous term plays a role mainly at the position of the 
shocks, where it is counterbalanced by the nonlinear term. The equality of 
these two terms leads to 

u = Au.5 (17) 

where Au is the velocity jump across the shock, and 5 is the shock thickness. 
There are thus two scales here: a large scale L corresponding to some average 
distance between shocks, and a dissipation scale 8 ~ u, very much smaller 
than L when v goes to zero. Distances away from both extremes make up 
the inertial range. 

In terms of multifractal language, the Burgers equation ( one averages, 
in the limit v — > 0, over an ensemble of initial states, or considers stochas- 
tic forcing on large scales) shows extreme multifractality, a situation called 
bifractality in the literature^]. The behavior of u is essentially linear be- 
tween shocks (u ~ x), and thus here h = 1,D(1) = 1. At the shocks 
themselves h = 0,D(0) = 0, since the shocks are discontinuities of the ve- 
locity field occurring at a point (in the v — > limit). The velocity variation 
across the shock is independent of distance, and the probability of being 
within a distance r is linear in r (cf. equations (12) and (13) for the case of 
1 dimension). 

There are thus two possible values for the exponent ph + 1 — D(h) (cf. sec- 
tion 1. 3.), namely p or 1, and therefore the dominant exponent £ p (equation 
(15)) characterizing the behavior of the structure functions in the inertial 
scale, is such that 

C P = 1, p>l (18) 

This is an extreme case of multifractality ( all exponents have the same 
value for integer p greater than 1), very much different from the case of 
three dimensional homogeneous, isotropic turbulence where the experimen- 
tally determined exponents remain relatively close to the scaling ones, which 
increase linearly with p (see equation (6)). 

However - as we have discovered - there is a whole range of multifractal 
behavior as the spatial extent of the stochastic force in the Burgers equation 
varies, and the situation is much more interesting. 
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II. 2. Stochastic forcing. 



For the stochastic forcing in (16) we take a Gaussian, such that in k 
space 

< f(k,t) >=0 

< f(k, t)f(k', t') >= 2D \kfS k> ^S(t - t') (19) 

The exponent f3 determines over which scales the forcing acts. For [3 > 
it acts effectively on small scales, whereas as (3 becomes negative, larger 
and larger scales matter. The limit relevant to forcing in three dimensional 
turbulence is that of large scales, of the order of the system size L. 

The range of values of (3 goes from /3 = 2, which corresponds to thermal 
noise, to (3 = —3/2. For values smaller than the latter, the statistics of the 
velocity field is independent of /3, unchanged from its behavior at (3 = —3/2. 
At (3 = —3/2 the system behaves as the steady state of the plain Burgers 
equation: it exhibits the extreme multifractal behavior discussed in ILL, 
characteristic of a shock dominated velocity field. For (3 > however, the 
presence of noise on small scales prevents the shocks from developing, and 
therefore the behavior appears chaotic, i.e. random and structureless. Thus 
as (3 moves from positive to large negative values, the velocity field goes from 



a chaotic to a shock dominated state, through an intermediate region|12] 
(—3/2 < /3 < 0), where for —1 < (3 < it displays complex dynamics of 
appearing, interacting and disappearing shocks. This region is one of rich 
multifractal behavior, and is the principal object of our study. It is through 
this region that one approaches the hydrodynamic limit of large scale forcing 
from a purely chaotic state. 

To be complete, we mention that for positive values of (3 one can use a 
renormalization group approach. As soon as (3 becomes negative, all sorts 
of non-linear terms become important in the equations, and the perturba- 
tive renormalization group approach breaks down. This approach has been 
usually applied |l|] to the equivalent KPZ (Kardar-Parisi-Zhang) equation 
for fluctuations of an interface height h(x, t), related to u by u = d x h. With 
a noise of the form considered, the renormalization group has also been 
applied to the Navier-Stokes equation |Q|. 

For (3 positive , close to zero, the scaling analysis leads to the following 
result for the exponents z and which appear in the scaling form assumed 
for 52 (r, r) =< (u(x + r,t + r) — u(x, t)) 2 >, namely 52 (r, r) = r^ 2 g(r/r z ): 

z + C 2 /2 = 1 (20) 
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and 



C 2 - z = -1 - p 



(21) 



The first relation is a consequence of Galilean invariance, the second of the 
fact that the coefficient Dq of noise fluctuations is not rescaled because of 
the non-analytic form of the noise. One obtains from (20) and (21) that 
( 2 = -2/3/3 and z = 1 + 0/3. 

We will from now on consider the region of negative /3, which is so to 
speak the gateway to hydrodynamic behavior. 

III. Three dimensional turbulence and the stochastic Burgers 
equation. 

We believe that because of the mathematical similarity of the Navier- 
Stokes equation with forcing, and the stochastic Burgers equation, the latter 
can be used as a key to the understanding of a number of issues in the 
statistical behavior of isotropic, homogeneous turbulence. In the work we 
have been doingQ, |j| ||, fj], 0], we highlight this similarity on a number of 
occasions, in different situations. To give a simple example here, we compare 
the Von Karman-Howarth relation for S3 for both equations. 

For the Navier-Stokes equation with forcing /, this relation takes the 
following form for the (equal time) 3rd order structure function S%j =< 
(u\ —U2) 2 (uij — U2j) >, where "1" refers to the point x + f, "2" to the point 
x, and "j" denotes the j-th component of u 

\d rj S Zj {r) = uAS 2 (r) - 2 < e > + < (ui - u 2 )-(fi - / 2 ) > (22) 

where S^r) =< {ui — U2) 2 >> while for the stochastic Burgers equation, 
where Ss(r) =< (u% — U2) 3 >, it reads 

\dS z (r)/dr = ud 2 S 2 /dr 2 - 2 < e > + < (m - « 2 )(/i - / 2 ) > (23) 
o 

The structural similarity of the two equations is clear. 

One can derive the above two Von Karman-Howarth relations in a straight- 
forward way from the space and time dependent S2, using the homogeneity 
in time of expectation values. More precisely, one writes that dS2(r, r) /dti + 
dS2(r, r)/dt2 = 0, where r = x\ — X2,t = t\ — This derivation highlights 
the fact, which we have several times pointed out in our work, that it is often 
useful for deriving equal time correlations to pass through time dependent 
calculations. Many identities can be obtained this way. 
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The two equations (22) and (23) are very similar. The 3 dimensional re- 
sult contains Kolmogorov's "4/5th" law for the longitudinal structure func- 
tion. In both cases < e > represents the energy dissipation rate. Since r 
belongs to the inertial scale the term multiplied by v is negligible in both 
equations in the zero viscosity limit. The noise dependent term can be 
evaluated in the equal time limit with the help of the Novikov-Donsker 
formalism |l5| . When the noise is cut-off at large scales (the hydrodynamic 
limit) this term leads to a subdominant correction of order (r/L) 2 . We will 
discuss later, for the stochastic Burgers equation, the general case when the 
noise ranges over small scales as well. 

Though this comparison of the Von Karman-Howarth relations is based 
on a simple case, we have found that the same similarity term by term, with 
an obvious display of the 3 dimensional space indices, holds for any other 
equation we have derived involving velocity or dissipation rate correlations, 
with the exclusion of course of terms involving pressure. 

We will discuss in the following three main points: 

(i) first, we are going to face for the stochastic Burgers equation the 
problem of turbulence, namely calculate, for small p, in the multifractal 
region (— 1 < (3 < 0) the exponents ( p characterizing the statistical behavior 
of velocity structure functions, 

(ii) second, we are going to give the general equation satisfied by the 
equal time correlation of the dissipation rate, and connect its intermittent 
behavior, which exhibits a hierarchy of exponents, to the intermittent be- 
havior of the velocity structure functions, 

(ii) third, we investigate the dynamics of the second order structure func- 
tion, and show how - even in the absence of any average flow - 5*2 satisfies a 
wave equation with characteristic velocity y/< u 2 >. These dynamic consid- 
erations enable us to disentangle, in our Eulerian framework, the intrinsic 
dynamical and the kinetic, ballistic characteristic times which describe the 
time evolution of flow structures. 

III.l. Multifractal exponents. 

We are interested in the region where — 1 < /3 < 0. Here also exists 
the possibility of scaling behavior, in the same way as there is Kolmogorov 
scaling for three dimensional turbulence, where the dimension of < e > or 
equivalently Dq, determines the dependence on distance of the S p 's in the 
inertial range. One thus has 

S p (r) ~ (D /Ly/ 3 r-^ 3 (24) 
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which corresponds to C, p = —p(3/2> and h = —{3/3. This is the value of h in 
the scaling regime. (Notice that at (3 = — 1 the exponents are the samejlC] 
as those of Kolmogorov scaling, equation (6).) 

This scaling regime is however dominant only in the region of (3 nega- 
tive close to zero, and gives way to multifractal behavior as (3 goes towards 
— 1. We are going to study this behavior directly on equations for the struc- 
ture functions derived from the stochastic Burgers equation. We proceed 
systematically discussing first 52 and S3, and then £4, £5 and general S p . 

(i) S 2 and S3. 

One cannot derive directly from the stochastic Burgers ( or from the 
Navier-Stokes equation in three dimensions) a closed equation for the equal 
time structure function S 2 . We therefore check numerically that 62 (r) be- 
haves in the following way 



S 2 (r) ~ (r/L)- 2 ^ 3 



(25) 



for all —3/2 < [3 < 0. Precise numerical results, and therefore a precise 
value of the exponent, can be obtained from evaluating the energy spectrum 
(E(k) ~ |£;|~ 1+2/3 / 3 ), related to 5*2 by Fourier transform, rather than from 
S 2 itself ( Figure 1). S 2 (r) thus scales, in the sense that (2 = —2/3/3 has its 
scaling value (cf. equation (24)). 



ho 




Figure 1: Graph of log E(k) as a function of logfe, where the energy spec- 
trum is E(k) ~ \k\~ 1+2/3 / 3 , for f3 = —0.8. The straight line for small k, 
drawn for comparison, has a slope of —1.53, which is the value of — 1 + 2/3/3 
at the given (3. 

As to Ss(r), it is determined from the Von Karman-Howarth relation, 
equation (23). In this equation the noise term takes in the equal time limit 
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(Novikov-Donsker formalism |[5|]) the form 

< («i - u 2 )(fi - f 2 ) >= 2(1/L 2 ) D \kf{l - coskr) (26) 

k 

The term proportional to "1" in (1 — coskr) cancels the — 2e in equation 
(23) because (1/L 2 ) J2k AdI^I' 3 is the total rate of energy input. One thus 
obtains from equation (23) (in the v — > limit) 

ldS 3 /dr = —2(1/L 2 ) V D \ kf coskr (27) 
6 , 

k 

The "coskr" term leads by rescaling to the following result 

S 3 (r) ~ r~P (28) 

for — 1 < (3 < 0, in the case where the noise does not have a cut-off at scales 
of order L. ( At j3 = —1 there is an additional logarithm, S3 ~ rlogr.) 

The exponents characterizing the inertial range behavior of S 2 and S3 
have therefore their scaling values throughout the domain — 1 < (3 < 0. For 
S 2 the result is based on simulations, for S3 the expression of the exponent 
is obtained from the Von Karman-Howarth relation. 

(ii) S4, S5 and general S p . 

For p > 4 scaling no longer holds through the entire — 1 < (3 < range. 
The following are the equations we obtain from the stochastic Burgers equa- 
tion after isolating the terms which in the inertial range go to zero when the 
viscosity does, and simplifying the noise terms 

\dS±{r)/dr = ^ud 2 S 3 /dr 2 - 2 < (e x + e 2 )(t*i - u a ) > (29) 



±dS 5 (r)/dr = l-vcPS^/dr 2 " ^2 E D \kfcos(kr) < ( Ul - m) 2 > 
~\[< (ei + e 2 )(m -u 2 ) 2 > 

- <{ei + e 2 )><{u 1 -u 2 ) 2 >} (30) 



dS p (r)/dr ~ ^2 ££> |Af coa(fcr) < («i - « 2 ) p_3 > 

A: 

+... < (ei + e 2 )(«i -u 2 ) p - 3 > (31) 
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The right-hand sides of equations (29) and (30) contain terms ( not written 
for equation (31)) which go to zero in the small viscosity limit, a noise 
dependent term and a dissipation rate dependent term. The noise term has 
the general form 

J2 D \kfcos(kr) < (m - u 2 r~ 3 >~ ^ls p - 3 (r) (32) 

k 

since dS^{r)/dr ~ J2k D^\k\^ coskr (cf. equation (27)) 

Therefore scaling behavior in S p is present, whether dominant or sub- 
dominant, whenever there is scaling behavior in Sp-%. Thus the presence of 
a scaling term in S2 , 53 and S4 guarantees the presence of one in any S p for 
p > 4. We have already pointed out that both S2 and S3 scale through the 
domain — 1 < (3 < 0. The case of S4 is trickier because of the absence of an 
explicit noise term in equation (29). We discuss it below. First we turn to 
extracting the multifractal behavior of S4 and higher order structure func- 
tions. This behavior becomes relevant when the associated exponents are 
smaller than the scaling ones, and therefore the corresponding non-scaling 
term dominates over the scaling one, since r/L <C 1. 

We first note that in /c-space both S3 and £4 depend on < u(ki)u(k2)u(ks) > 
,ki + k2 + k% = 0, the first one through its definition, the second one through 
the e dependent term in (29). We thus make the following general ansatz 

|&i| m |£:2Mfc3| M3 

Im < u{ki)u(k 2 )u{k 3 ) >~ — - — h permutations (33) 

K1K2K3 

The constraint that Ss(r) ~ r~ /3 (cf. equation (28)) leads to u\ + U2 + Us = 
1 + (3. We can show that the lowest exponent is obtained when u\ = u 2 = 
U3 = u/3 = (1 + /?)/3. Putting the ansatz into the 2nd term of (29) leads to 

dS/±/dr ~ v f dadk\dk 2 dk 3 sin(fc 1 r) ■ eX p —ia(k\ + k2 + ks) 

(34) 

Performing the k integrals with a cutoff 5 and then integrating over a, with 
< ji\ < 1, one obtains 

dS A /dr ~ u(2TT/5f 2+ ^(l/5)r-^ (35) 

and thus, with ^i = ^2 = 

S *(r) - sm^r^ 3 (36) 
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It is important to note here that the non-scaling behavior arises from 
the term in the equation which involves e. The expression for S4 contains 
two results: 

(i) the fact that in the limit v — > , 

v „ J l+2(l+/3)/3 (37) 

whereas in the scaling limit v ~ J 1- ' 3 / 3 . (By writing that at the dissipation 
scale 5, the characteristic eddy time tg ~ 5/5 h is of order of the dissipation 

time 5 2 /v, one finds v ~ 5 1+h ) One thus has a new dissipation scale in S4, 

1 

namely 5 ~ v Y+h ± . This dissipation scale depends on the corresponding mul- 
tifractal exponent /14 = 2(1 + /3)/3. For the dominant term this multifractal 
exponent has to be construed as the one which minimizes £ p (cf. (15)). 
(ii) second it gives the non-scaling exponent £4 = (2 — /3)/3, which being 
smaller than the scaling exponent £4 = —4/3/3 in the region —1 < /3 < —2/3, 
dominates over the scaling term in this region. 

We now have to get back to the question how scaling behavior arises 
in S4. One can show that it arises through the vdSil&r contribution in 
S3 present in the Von Karman-Howarth relation (cf. equation (23)). It 
corresponds to fi± + ^2 + ^3 = 2 + 2/3/3 in the ansatz for S3 (see above) with 
however / \ii = fJ>3- 

One can now proceed along the same lines to find the behavior of Ss(r), 
taking as a starting point an ansatz similar to the one used for S4, but now 
for < u(ki)u(k2)u(k^)u(k4 : ) >,k\ + k% + k% + foj = 0. There are now four 
//'s, the sum of which is constrained by the known behavior of S4 in two 
different regions —1 < j3 < —2/3 and —2/3 < (3 < 0. We know already that 
in S5 because of the presence in equation (30) of the noise term, a scaling 
contribution will be present. The question that is to be settled through 
making the ansatz on the 4-point function, is whether there are regions in 
which the scaling term is subdominant, as happens for S4. The answer is 
yes, and one finds that there are three different regions: 

(i) —1/2 < /3 < 0, where scaling behavior dominates, and thus £5 = —5/3/3, 

(ii) —2/3 < P < —1/2, where S5 does not scale, £5 = (3 — 4/3)/6, and this 
exponent is smaller than the scaling one and therefore the corresponding 
term dominates in Ss(r), 

(iii) — 1 < P < —2/3, where S5 has still another multifractal exponent, 
£5 = (5 — /3)/6, which gives the dominant behavior in this region of /?. The 
three exponents connect smoothly at the end points of each interval. In each 
interval all three terms are present, but the term with the smallest exponent 
dominates. The first four £ p 's are shown || in Figure 2. 
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-2 -1.5 -1 -0.5 

P 

Figure 2: Exponents (2,(3X4 an d Cs vs - P, for — 2 < (3 < 0. The lowest 
curve is that of (2, the others above it are in order (3, £4, £5. 

The following general scenario emerges from these results: as p increases, 
simple scaling with ( p = —p(3/2> occurs over a progressively diminishing 
range of values for [3 close to zero (and negative). Over most of the consid- 
ered domain therefore, multiscaling occurs as soon as p > 4, with the ( p 's 
continuous and piecewise linear, the number of linear segments increasing 
as p gets larger. As (3 — > — 1 all the £ p 's for p > 3 go towards 1. This 
extreme multifractal regime is a manifestation of the increasingly important 
role played by shocks as the noise acts on larger and larger scales. 

Several remarks are in order here: 

(i) if one extracts a fractal scaling exponent for velocity variations from 
the calculations, as we have done above for S4 (equations (12) and (37)), one 
finds a different value for /15 in each of the three regions of /3, where different 
Cs's dominate, namely h 5 = -(3/3 for -1/2 < (3 < 0, h 5 = 1/2 + 2/3/3 for 
-2/3 < (3 < -1/2, and h 5 = (1 + (3)/Q for -1 < < -2/3. Thus h 5 is 
continuous and piecewise linear, and goes to zero as (3 — * —1, which is a 
reflection of the increasing dominance of shocks. The same is true for all 
Zip's with p > 4. 

(ii) one can also calculate continuous and piecewise linear fractal dimen- 
sions D(h p ) with the help of equation (15), assuming that the corresponding 
hp minimizes the right hand side, and using the values of h p and ( p which 
result from the " ansat^ calculation. One finds that all fractal dimensions 
tend towards zero as /3 — ► — 1, which again is consistent with the dominance 
of shock structure. 

(ii) we cannot show in general that our calculation based on an ansatz 
in /c-space, and the assumption of the equality of /i's in S5 (cf. equations 
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(33) and (34)) leads to the "true" dominant behavior in each domain. It is 
possible that our continuous, piecewise linear ( p 's, are only an approximation 
to the "true" function Cp(P)- 

III. 2. Dissipation rate correlation and intermittency. 



By studying the full equation satisfied by the dissipation rate correlation 

< e(x + r)e(x) >~< e > 2 {r/L)~^ (38) 

we are able to find expressions for the intermittency exponent u in terms 
of static and dynamic exponents of velocity field correlations. Here e(x) = 
u(du/dx) 2 for the Burgers equation and e(x) = %(diUj + djUi) 2 for the 
Navier-Stokes equation are the local dissipation rates. In our previous dis- 
cussion, we have taken the energy dissipation rate e to be a constant, and 
this is all that is required to obtain Kolmogorov scaling of the structure 
functions. In this section e(x) is considered to be a fluctuating quantity 
which has non trivial correlations, as experiment shows. One still has 
< e(x) >= e = constant because of homogeneity. 

The following two relations have been proposed for the intermittency 
exponent \i\ 

Mi = 2 - Ce (39) 

and0, [17] 

M2 = 2C 2 -C4 (40) 

The first one, the most discussed, because experimentally the value of ^6 ra 
1.8 agrees with that of u « 0.25|l8[], is essentially obtained by a scaling 
argument, which uses the dimension of e, namely V 3 /L, to set < e(x + 
r)e(x) >~ S G {r)/r 2 ~ (r/L)^. 

The advantage of our approach lies in the fact that relations between u 
and structure function exponents £ p are derived directly, and simultaneously, 
from the equation satisfied by the dissipation rate correlation. This equation 
can be derived from the stochastic Burgers or the Navier-Stokes equation by 
considering correlations in both space r and time r, and then passing to the 
r — » limit. In this limit the noise term can be expressed using the Novikov- 



Donsker formalism] 15]. One finds in this way, with e\ = e(x + r, t + t), €2 = 

e(x,t) 

< eie 2 > = \d T < (ei - e 2 )(iti - u 2 ) 2 > -\d r < (ei + e 2 )(wi - u 2 f > 
4 b 

--d r < (ui - u 2 ) 2 {e 2 U2 - eiui) > +^<9 2 < (ei + e 2 )(u 1 - u 2 ) 2 > 
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+^ < («i - «2)(e 2 /i - ei/a) > (41) 

The 3rd term on the right-hand side ensures Galilean invariance together 
with the first term (the left-hand side is Galilean invariant). The viscosity 
dependent term, which is connected to d 3 S^/dr 3 (cf. equation (30)), goes 
to zero for inertial r in the zero viscosity limit. 

In order to show again the mathematical similarity of expressions de- 
rived from the Burgers and Navier-Stokes equations, we show the equivalent 
expression in three dimensions derived from equation (2): 

< ei£2 > = ^<9 r < (ei - e 2 )(ui - u 2 ) 2 > -jd Tj < (ei + e 2 )(itij - u 2j )(ui - u 2 ) 2 > 

+ \ dr o < ("! ~~ ^) 2 (eiuij - e 2 u 2j )) 

< (ei + e 2 )(ui - u 2 ) 2 > +-d n d rj < eiu 2i u 2j + e 2 uiiUij > 

+^ < (uu - u 2i )(e 2 f u - eif 2i ) > (42) 

Apart from the pressure term and a more complicated viscosity term due to 
the difference in structure of the definitions of e in the Burgers and Navier- 
Stokes case (see the beginning of this section), the two equations correspond 
to each other term by term, with an obvious generalization of space indices 
when going from one to three dimensions. 

Now going back to equation (31) withp = 6, one sees that the expression 

< (ei + £2X^1 — u 2 ) 3 >, which occurs in (41), is precisely the term in dS^/dr 
which, as argued in section III.l., leads to intermittency. Therefore from 
(41), < t\€ 2 > (in the r — > limit) contains the intermittent behavior 
(r/L)-^ 1 , with 

Mi = 2 - Ce (43) 

as given in equation (39). 

As to the first term on the right hand side of (41), one can showQ that 
the expression < (ei — e 2 )(u± — u 2 ) 2 > appears in dS^/dr, where it is the 
only one involving the dissipation rate, and therefore leads to intermittency. 
There is thus a contribution here to the intermittent behavior of < £162 > 
of exponent 

m = z A , 2 - C4 (44) 

where z^ 2 characterizes the behavior of the second order partial derivative 
of S4 in time, in the limit r — » 0. The origin of fj> 2 is thus dynamical. If 
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simple scaling in time holds, then 2:4 2 = 2z, where z = 1 — h, with the value 
of h equal to its scaling value, z here is the dynamical exponent, not the 
"frozen turbulence" exponent of value 1, which characterizes the advection 
of small structures by large ones. Our preceding result and remarks apply 
as well to Navier-Stokes turbulence. In this latter case z = 2/3, which is 
numerically equal to C2 (we are going to show in III. 3. that this result is 
general and exact). Substituting £2 for z (recall that in the scaling limit 
-24,2 = 2z) in (44) leads to the result given in equation (40), which thus 
appears as a static approximation to what our derivation shows to be the 
dynamical intermittency exponent given by equation (44). 

For the Burgers equation the two intermittency exponents of equations 
(43) and (44) are the two main ones. For the Navier-Stokes equation we can 
only assert that these same two occur as well, because our discussion does 
not take into account the pressure term in equation (42). 

III. 3. Dynamic behavior. 

Except for the discussion of ^2 m the preceding section, our concern up 
to now has been with the equal time structure functions. We now address 
the problem of their dynamical behavior. We are interested in relationships 
between dynamic and static exponents, and also in shedding light on Taylor's 
frozen turbulence hypothesis in the case when there is no average flow field. 
In particular we wish to understand how it happens that the square root of 
the rms fluctuations of the velocity field replaces the average velocity when 
the latter is zero, thus allowing ballistic behavior with z = 1 (z is defined 
by t ~ r z ). The objects of our study are now the space and time dependent 
structure functions 

S p (r,T)=<( Ul -u 2 r > (45) 

where u\ = u(x + r, t + r),«2 = u(x,t). The generalization to the three 
dimensional case is straightforward. 

We will concentrate on S2 ■ One can derive the following equation from 
the stochastic Burgers equation Q 

dS 2 (r, t)/8t = ^8T 3 /dr+ < u x f 2 >-< u 2 /i > (46) 

where 

T 3 (r,r) = - < {u l + u 2 ){u l - u 2 f > (47) 

which apart from additive constants is the same as < u\u2 + >• The 
term on the left-hand side and the first term on the right-hand side form a 
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Galilean invariant pair. In the r — ► limit T3 does not contribute because 
of symmetry reasons. In this limit there is a discontinuity in the noise term 
because < 1x1/2 > contributes for r > 0, and < U2/1 > f° r t < 0. One thus 
has, using equations (23) and (27), 

dS 2 (r,r = 0+)/cV = (1/L 2 ) D \kf coskr = -—dS 3 /dr (48) 

k 12 

Assuming simple dynamic scaling for the first time derivative of S2 (in the 
r — > limit), with r ~ r z , equation (48) leads to the following relation 

C2 - z = C3 - 1 (49) 

This equation is the same as equation (21). However here it follows from 
an exact equation, whereas before it was obtained from a renormalisation 
analysis. Moreover z here is precisely defined as the exponent which char- 
acterizes the behavior of the first order partial derivative of S2 in time in 
the limit r —> 0. 

Since £3 is known from the Von Karman-Howarth relation (equations 
(22) or (23)), this equation relates the temporal and spatial exponents which 
characterize the behavior of the 2nd order velocity structure function. Since 
C3 has its scaling value set by the Von Karman-Howarth relation, any scaling 
violations in £2 has to be compensated by an equal one in z. Introducing 
the value of (3, one thus has in the case of the Burgers equation 

C2 - z = -(3 - 1 (50) 

and in the case of Navier-Stokes 

C2 - z = (51) 

Thus £2 and z are not independent, the knowledge of one determines the 
other. This is the first constraint we have found for £2, for which none can 
be found when one limits one's investigations to static quantities only. In 
particular, in the Navier-Stokes case C2 = 2/3 = z, whereas in the Burgers 
case one obtains z = 1 + (3/3. The latter results are consistent with the 
simple Kolmogorov type scaling argument which entails z = 1 — h. 

As r — > what matters is clearly this dynamical z, the one appropriate 
for a Galilean invariant situation. However as soon as r departs from zero, 
the ballistic behavior with z = 1 asserts itself. We have checked this numer- 
ically for S2(r = 0, r) and S^{r = 0, r), for which, if dynamical scaling holds 
and for example S^r, r) = r&g(r/r z ), time dependence is of the form t^ 2 ^ z , 
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Figure 3: Log%oSi(r = 0,r) vs. logiQT for a) (3 = —0.5 with a dashed line 
of slope C<i/ z = 2/3 with z = 1, and b) /3 = —1 where the dashed line has 
a slope of 0.92 close to the numerically observed value of (4. The expected 
slope is Ca/z, and the numerical results allows one to distinguish between 
z = 1 and z = 2/3, the value of z = 1 + /?/3 for (3 = —1. 



and similarly for S4. Numerically one is able to distinguish satisfactorily 
between the dynamic and ballistic values of z (Figure 3). One thus verifies 
that as soon as r is positive, ballistic behavior with z = 1 occurs. 

The question now arises in which way ballistic behavior emerges, and 
with it the use of Taylor's frozen turbulence hypothesis, in the case when 
there is no average flow, i.e. < u(x, t) >= 0. 

In reference [Q] we have shown that if one differentiates relative to r 
equation (46), one is lead to the following equation 

d 2 S 2 {r,r)/dT 2 =< u 2 > d 2 S 2 /dr 2 + (52) 
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The term on the right-hand side is a result of the fact that 



dT 3 /dr oc< u 2 > dS 2 /dr (53) 

after use of the assumption that in the v —* limit the term < {u\ — 
U2) 2 {u\ 2 + U2 2 ) >~ 2 < u 2 > $2- The latter assumption arises from the 
observation already made by Tennekes[19] that large scales eddies advect 



inertial scale information past an Eulerian observer. Here we show that 
this assumption is encapsulated in the fact that 52 (r, r) satisfies precisely a 
wave type equation of characteristic velocity given by the rms fluctuations 
of the velocity field. One expects this behavior to occur over time scales 
large compared to the dissipation time and small compared to the turnover 
time of the large scale structures in the system. A detailed discussion of the 
other terms occurring in the equation can be found in reference [pi] . 



IV. Remarks on intermittency. 

Before embarking on these remarks one should point out that the nature 
of turbulence is different for the Burgers and Navier-Stokes equations: for 
example vortex stretching is believed to be an important ingredient in three 
dimensional developed turbulence. 

Intermittency - the non-scaling behavior of the structure functions in 
the inertial range - is a hallmark of three dimensional turbulence. The 
language of multifractality is a convenient way to describe it. What is the 
origin of intermittency in the statistical behavior of turbulence? The answer 
is not clear, though intermittency has been connected to the presence of 
vortex filaments in the flow. In one experiment fl20|] , where the size of the 
filaments is several times the dissipation scale, they are associated with 
events in the velocity field where the velocity derivative has large jumps. 
This is of course what happens across shocks, which play the role of coherent 
structures in the one dimensional stochastic Burgers equation. Here one has 
a clear connection between intermittency and the presence of shocks, though 
we are unable to give a numerical measure of the number and sizes of shocks. 
Typically the velocity variation across a shock occurs on length scales of the 
order of the dissipation scale. For (3 negative close to zero, shocks are barely 
apparent in the velocity profile, and the structure functions show scaling 
behavior. As f3 approaches —1 the shocks play a larger and larger role, and 
intermittency, the difference between the actual values of the £ p 's and their 
scaling values, increases correspondingly (for p > 4). At j3 < —3/2 the 
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shocks are present in full, dominating the velocity profile, and intermittency 
is extreme: all £ p 's are equal to 1. There is thus an obvious link between the 
dynamics of shocks - the small scale coherent structures - and intermittency. 
We provide two other insights: 

- we connect - not by a self-similarity argument, but from the exact 
equation - the values of the exponents measuring intermittency in the energy 
dissipation rate to those measuring intermittency in the velocity structure 
functions (see III. 2.), 

- we show that in the equations for the velocity structure functions the 
terms responsible for intermittent behavior are those which contain the en- 
ergy dissipation rate. Intermittent behavior at the inertial scale is thus a 
consequence of dynamics which occurs at dissipation scales (see III.l.). 

For the stochastic Burgers equation we are of course able to provide an 
extra bonus: namely, with the help of an ansatz, we are able to calculate from 
the basic equations the low order structure function exponents as (3 varies. 
Such a calculation remains the "holy grail" for statistical three dimensional 
turbulence. [21] 
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